A dynamical mean-field theory approach 
to superconductivity and 
antiferromagnetism in a strongly 
correlated electron 
system. 



H. Watanabe^ and S. Doniach^'^ 
Departments of Applied Physics^ and Physics^ 
Stanford University, Stanford, CA 94305, U.S.A. 

February 1, 2008 

Abstract 

In this paper, we present the results of numerical studies of superconduc- 
tivity and antiferromagnetism in a strongly correlated electron system. To do 
this we construct a Hubbard model on a lattice of self-consistently embedded 
multi-site clusters (in practice two sites) by means of a dynamical mean-field 
theory in which intra-cluster dynamics is treated essentially exactly. We show 
that a class of characteristic features which have been seen in the excitation 
spectra of high-Tc cuprates including the pseudogap and the spin-flip resonance 
feature seen in neutron scattering studies, as well as their interplay with the 
onset of a pairing correlations, can be captured within a dynamical mean-field 
theory in which short-wavelength dynamics are rigorously treated. Thus we 
infer that the observation of the neutron scattering resonance in the supercon- 
ducting state of the cuprate superconductors does not appear to be directly 
tied to their quasi-2D character. 

Although our approach is defined strictly in terms of fermion degrees of free- 
dom, we show that we can readily identify the emergence of effective low energy 
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bosonic degrees of freedom in the presence of a well-defined broken symmetry 
phase as long as their dynamics are dominated by short-range, short- wavelength 
fluctuations. Our exact calculations reveal that the dynamics of the spin de- 
grees of freedom and the onset of superconductivity are strongly entangled. In 
particular, the dynamics of staggered spin degrees of freedom builds up coher- 
ence and a resonance-like sharp feature as superconducting pairing correlations 
set in (this feature diminishes in the normal phase). At the same time a spin 
gap develops in the staggered spin susceptibility. Under conditions of super- 
conducting broken symmetry our approach thus extends static BCS mean field 
theory to provide an exact treatment of quantum fluctuations of the BCS or- 
der parameter within self-consistent dynamical mean-fleld theory. We find that 
both equilibrium and dynamical properties of our model can provide a consis- 
tent interpretation of experimental observations. 

I. INTRODUCTION 

Understanding the physics of strongly correlated electron systems has been one of 
the most difficult challenges in condensed matter physics for several decades. The im- 
portance of electron-electron correlations in high-Tc superconductors was recognized 
soon after its discovery and much progress has been made in understanding its 
phenomenology. As a result of intensive efforts in both theoretical and experimental 
studies, the detailed nature of the microscopic mechanism for high-Tc is starting to 
become clearer and the idea that the superconducting instability is driven by Coulomb 
interactions is becoming more widely accepted, although it still remains controversial. 

Besides undergoing a superconducting transition at high temperatures, the cuprate 
superconductors exhibit other features characteristic of strongly correlated systems 
including a pseudogap, spin fluctuation resonance, and non-Fermi liquid behavior in 
the normal state, some of which are also found in other transition metal oxides. 

The key motivation of the present paper is to use an extended local dynamical 
mean field approach to study the effects of strong correlations on the superconducting 
instability. By going beyond the static mean field BCS treatment we are able to show 
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that coulomb correlations can stabilize the superconducting state and also lead to 
a relationship between the onset of superconductivity (or more generally, pairing 
fluctuations) and the spin dynamics similar to that seen experimentally. We are 
also able to study the onset of the pseudogap and its effects on the single particle 
properties of the system. 

Considerable experimental evidence has accumulated which suggests that the mi- 
croscopic pairing mechanism of high-Tc superconductivity may be already manifested 
at a relatively short length scale. Both superconducting coherence and spin-spin corre- 
lation lengths in high-Tc cuprates are relatively short. The in-plane superconducting 
coherence length, for example, is estimated to be roughly ~ 15A which indicates that 
Cooper pairs, on average, span only a few lattice spacings. Also, the correlation length 
for spin degrees of freedom estimated from a relatively broad resonance peak width in 
a momentum space is also of the order of a few lattice spacings 0, ^ H, ^, |[. In ad- 
dition, STM measurements of the influence of magnetic impurities on local electronic 
structure of high-Tc cuprates clearly suggest that Cooper pairs are local entities, and 
can exist in a microscopically confined region ||^. 

These measurements suggest that a theoretical description based on a real space 
representation(in contrast to the BCS momentum-space representation) in which fluc- 
tuations in a relatively localized region are rigorously treated can be a reasonable 
starting point to study the superconducting instability. 

In this paper we show that a cluster-based dynamical mean-field theory approach 
is particularly useful for this purpose. Based on this approach, we are able to study 
the nature of short-ranged quantum fluctuations of a strongly correlated system in the 
presence of well-defined superconducting and antiferromagnetic broken symmetries. 

The numerical results we report in this paper are an extension of earlier work|P] 
(unpublished) and are based on a Hubbard-like lattice model consisting of two-site 
clusters using the self-consistent dynamical mean-field theory approach of Kotliar 
and Georges^]. In the dynamical mean-field theory approach, a lattice problem (for 
which an exact solution is generally difficult to obtain) is mapped onto an impu- 
rity problem coupled to an effective bath (which is typically much more tractable) 
- thus the advantage of this theory is that once the model is constructed it can be 
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solved in principle exactly without further approximations when supplemented by a 
self-consistency condition which is derived by requiring that the mean-field theory 
becomes exact in — oo limit (d: dimension) |jlO| . The crucial point here is that the 
effective bath is allowed to be time-dependent and thus highly non-trivial quantum 
nature of (local) dynamics is retained. This approach has offered new insights to 
the physics of strongly correlated electron systems which are difficult to obtain from 
pre-existing theories |n[] . 



Here, in order to address the quantum aspects of non-local but short-ranged 
fluctuations, we extend the original formalism to a self-consistent two-site cluster in 
order to be able to explore the superconducting part of the phase diagram. In order 
to treat superconductivity , we explicitly allow a U{1) gauge symmetry breaking for 
both diagonal (on-site) |jl2| and off-diagonal(nearest-neighbor) pairing channels 



inside a cluster. Within this model we are able to confirm the presence of a super- 
conducting phase for reasonable Hubbard-type parameters. By virtue of enforcing 
self-consistency, our model effectively mimics an infinite system, and thus is able to 
sustain a generic symmetry breaking. Thus, our model facilitates a study of short- 
range fluctuations in a well-defined broken symmetry phase such as superconducting 
or Neel order. 

Within this generalized model we are able to show that the onset of super- 
conductivity is intimately coupled to antiferromagnetic spin fluctuations and that 
a resonance-like feature is indeed seen to emerge in the dynamical spin susceptibility 
as the superconducting order parameter switches on. Because we are able to solve the 
local cluster problem essentially exactly, the presence of the Hubbard on-site repulsion 
U will automatically favor an off-diagonal (nearest neighbor) superconducting order 
parameter over an on-site one. We interpret this by analogy with BCS mean field 
theories of the cuprate superconductors 0, [l^ as favoring a d-wave like as opposed 



to an s-wave like order parameter symmetry. Recently Lichtenstein and Katsnelson 



T3[] and Maier et al [0 have developed approaches to dynamical cluster calculations 
for clusters with 4 or more sites which show properties similar to those found in our 
coupled 2-site model. Our calculations allow us to examine the spectra of excita- 
tions in our representation of the strongly correlated system which tie in well with 



experimental observations. The approaches of Lichtenstein and Katsnelson [13] and 



of Maier et al|]T6| may be expected to yield similar results. Recently, several groups 
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T7|, |T8| have argued that the resonance peak seen in inelastic neutron scattering may 



be thought of as a type of "spin - flip exciton". Our results suggest that the reso- 
nance is a generic property of the strongly coupled superconducting state and is not 
dependent on details of the cuprate band structure. 



The following is the organization of the paper. Section II will describe the details 
of our methodology and technical aspects. In section III, we present our results as 
well as the interpretation. In section IV, we discuss the overall aspects of the results 
obtained in section III and also some issues which were not addressed in the preceding 
sections. We then suggest some future work and close the section with a summary. 
Supplemental derivations and discussions are given in appendices. 



II. Methodology 



A. Model and Formalism 



A two-site cluster model is the minimal model needed to address the nonlo- 
cal nature of systems with order parameters for antiferromagnetism (staggered mag- 
netization) and superconductivity (nearest-neighbor pairing) and we will base our 
analysis on this model (Fig.l). We retain on-site Coulomb repulsion U, intra-cluster 
hopping t' and inter-cluster hopping t. (The reason we distinguish inter-cluster from 
intra-cluster hopping will become clear below.) Thus, the physics which our model 
represents is described by the following partition function: 



Z= Tr [e-f^^] 



<i,j>,a 
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i j 

+ V ^aIJe]^^ (1) 

<i,j>,(j 

+ {h.c.) (2) 



where is a bare energy level of the d-orbital of our model (i.e., Cu — Sd^^-y^ for 
high-Tc cuprates) and // is a chemical potential, d^l^ is the creation operator of an 
electron with spin a on the ith site of the A sublattice. < i,j > and repre- 
sent the sum over the inter and intra-cluster nearest neighbors, respectively, hz is 
an infinitesimal staggered magnetic field. 1] is an infinitesimal U(l) gauge symmetry 
breaking field for off-diagonal pairing, which we take to be real (and all the anomalous 
components as well). 



Now, we self-consistently embed our cluster into a lattice. Although the details 
of energetics of charge and spin degrees of freedom may depend on the lattice, most 
of the qualitative physics which appear to be manifested in real systems are shared 
by the tight-binding Hubbard-like lattice models embedded on a lattice which has 
a bipartite nature and a relatively smooth, structureless DOS [|Ty], In this 

study, we chose to work with a Bethe lattice. 

The self-consistency condition basically arises as a result of seeking a homoge- 
neous solution in the lattice problem. Here, we essentially follow the prescription 
developed by Kotliar and Georges [0, and extend it to the case of two-site 
clusters in the presence of superconductivity . The basic idea is to systematically 
expand (1) with respect to the inter-cluster hopping and pairing amplitude (scaled 
as t = r] = -^) and integrate out over the ligand degrees of freedom. Due to the 
above scaling, to order 0(1), only the lowest order term(the two-point propagator of 
the ligand multiplied by the probing fields) survives the integration and we obtain 
the following mean-field partition function 



DtfjDtlj^e-^'" (3) 

local 
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where 



S, 



eff 



[ [ dTdT'i^UT)G-\T -t'U(t') 

Jo Jo 



\){nBi{T) - ^; 



■0^ = {dAp dsp (i^j,! dBi) is a four-component Nambu-Gorkov spinor representing the 
local degrees of freedom. Go is a 4 x 4 cavity propagator for the cluster which can 
be thought of as a local noninteracting propagator of a particular cluster where only 
the on-site U for this cluster is turned off but all the other U are still present. From 
the above effective local action , a total propagator G{t) can be calculated as: 



G{t - r') 



■IMF 



(4) 



Once a total propagator is obtained, the cavity propagator in the Bethe lattice case is 
then given through the following self-consistency condition which can be conveniently 
expressed in frequency space: 



( T^^AAi'i'^n) T^ABii^n) ^AA(^^n) ^^Asi'i'^n) \ 

I^BAi'i^n) V]^Bi'i^n) A^j^Ai'^^n) 

AAAi'i'^n) A%{iuJn) V^J-^{iuJn) V^J-^iiuJn) 
V A^^Ai^^n) A^^Bi^^n) V^BAi^^n) V^BBi^^n) J 



where 

V']^^{iun) = iujn- ed + fi-^ - fG^^Bi'i^n) , 

AAAii^n) = t^Fl^B^lUJn), 



7 




and similarly for the other elements of G^^. Full expressions are given in the Ap- 
pendix. 

A few comments are in order. As a result of the decomposition of the lattice into 
clusters and treating the physics inside and outside the clusters in a different fashion 



hopping matrix elements no longer reflects the actual physics so that they should not 
be compared directly. Here, our main focus is on the short-range superconducting 
and antiferromagnetic spin fluctuations; i.e. that at part of the physics which intrin- 
sically depends on the short-range dynamics, should not depend qualitatively on the 
construction of the lattice, and should still be manifested in the intra-cluster dynam- 
ics. Considering the overall qualitative consistency of our results with experiments 
(please see Sec. Ill), we believe that the physics presented here is not an artifact due 
to the specific construction of our model. 

As will be shown below (see results) the self consistent numerical solution of 
the above model leads to a phase diagram exhibiting both antiferromagnetism and 
superconductivity. However, we find that the superconducting critical temperature 
is relatively low for us to be able to obtain reliable numerical results for the exci- 
tation spectrum of the model. In order to strengthen the tendency of the model 
to go superconducting , we perform a heuristic extension of the cluster expansion 
of the lattice action to higher orders in ^ and incorporate the effect of inter-cluster 
particle-particle correlations. For example, at O(^), by employing the idea of the 
Hartree-Fock theory, we can extract a inter-cluster particle-particle correlation from 
the following four-point propagators in a local effective action Seff'- 

rP rP rP rP . „ . 

/ / / / dridr2rfr3dr4doAT(n)'''(^oBi(T2)"^4Ai(T^3)t^oi?T(^4)G]j^AlBA(n,'r2;T3,r4) 

JO JO Jo Jo 



(see, however, ^), the ratio of the bare values of the inter- and intra-cluster 
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+ |F(0)| / / dTidT2doBi {Ti)doAi{r2)dA^ (ri)(isT(T2) 
Jo Jo 



where |-^(0)| is the average of equal-time inter-cluster anomalous Green's function. 
doAa and dAj represent local and neighboring cluster degrees of freedom respectively. 
In the above, since each local degree of freedom carries a pre-factor and we sum 
it over all d nearest neighbors, the above term becomes the order of O(^). 

We then incorporate the intcr-cluster pairing correlation from higher orders in ^ 
by employing the idea of the BCS theory. In order to preserve the effects of transla- 
tional invariance within the Bethe lattice model, we then substitute the inter-cluster 
pairing by the intra-cluster pairing correlations which can be rigorously calculated : 



_ginter ^ doAiadBj-adBj-adoAia 

— —g^"'^'" ^ < doAladB]-a > dsj-crdoAia + 9^"'^'^'^ X! doAi^dBj_„ < dBj-adoAia > 

- E < dlAljosl-a > dsj-jAia + 5'"*"' E ^^L^]-^ < dlsj-adoAia > 

= -fi^"'*" I ^AB^dBj-adAia + dA^^dB^.^isY 



where < doAi^doB^-a >= ^ab is the thermal average of instantaneous intra-cluster 
pairing which is calculated from the local quantum dynamics. For simplicity inter- 
cluster and intra-cluster coupling constants arc taken to be the same, g"^^''"- = g^^t^r 
and wc set their value equal to t. Wc then take a mean-field limit with respect to the 
intcr-cluster hopping matrix elements as well, i.e., C = and obtain the following 
equations: 
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V 



TT 
AA 



BB 



'^^Cab^ba'^bb I 



4.2 cU fit 
<A 

- t\e^j,F'^'B{^^^ 



it 

BA^ BB 



'j-^n)), 



(5) 



+ 
etc 



.2tTJ- cU pit I 



(6) 



where only representative components are shown and the full expressions are given 
in the Appendix. 



We will see that the effects of the additional inter-cluster coupling are manifested 
in all the components of the propagator. As we show in section III, superconductivity 
already appears in the absence of inter-cluster pairing correlation terms. When the 
additional inter-cluster coupling is included we find as expected that the inter-cluster 
coupling stabilizes the superconducting phase. In particular, Tc is enhanced for all 
the doping levels and the single-particle gap (i.e., the superconducting gap) becomes 
more pronounced. This enhancement in turn enables us to perform computations 
at higher temperature and hence to obtain more detailed excitation spectra in the 
superconducting region of the phase diagram. 



B. Numerics 

The local dynamics (3) and (4) is basically that of a degenerate impurity An- 
derson model(with the U(l) gauge symmetry breaking) and can be solved in various 
ways. Here, we chose to solve it numerically using the quantum monte carlo(QMC) 
algorithm of Hirsch and Fye Since there are two vertices in the local action 



(3), two Ising variables are introduced in our case. Then importance sampling is 
performed sequentially based on the ratio of the statistical weights for flipping Ising 
variables, which are given by 
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where (T^(B)(Tj) is the Ising variable at a time step rj(i=l,...,L; L:the number of time 
shces) which is introduced to decouple the Hubbard U term of A(B) site in a cluster. 
Here, we would like to warn the reader not to confuse Ising variables with electron 
spins. Ising variables are basically the fluctuating fictitious magnetic field which cou- 
ple to electron spins. A is a constant which is given by the relation, cosh{X) — e^. 
^AA('^i) "^j) is a diagonal component (in terms of imaginary time index) of the instan- 
taneous Matsubara Green's function at time step Ti of an arbitrary sweep, and should 
not be confused with the averaged Matsubara Green's function defined in (4). If 
the flip at the time shce is accepted, each component of the Matsubara Green's 
functions is updated according to the following rule: 



V 1 + (1 - G?^'^(rfc,r,))(e^^-(-^) - 1) ^ 



(7) 



^Vab{n) = -2\aA{Ti)\A]><A] 
-2\(jA{.n)\Ai><Ai 
-2X(jB{Ti)\B ]>< B ] 
-2\aB{Ti)\B i>< B i 
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a,b,c = {At Bt Al, Bl} 



where AVab{Ti) is a diagonal tensor which represents the amount of change in the 
action associated with the flip of an Ising variable at time slice r,. 

The above set of equations (3)-(5) is iterated until the convergence is achieved. 
In practice, a particular form of starting propagator is guessed depending on the type 
of the anticipated solution. The impurity problem (3) and (4) might be solved more 



efficiently by combining exact diagonalization and a Pade approximation [|TT]. The 
energy resolution of this approach basically depends on the number of energy levels 
to be used as a basis. Thus, it may not be the best approach when the spectral weight 
is widely spread and a tiny feature is under search (such as a dip structure in the 
ARPES spectrum in this case. Please see Section 8.2). Probably, this approach is 
best used when the shape of the spectral function is qualitatively known or reasonably 
guessed. 

A few remarks are in order regarding the sampling algorithm. A complete algo- 
rithm would be to integrate over the entire phase space of aA^Ti) and aB{Tj){i,j = 
1,2,..., La)- The size of the entire phase is 2^x2^ and it will be impractical to sample 
all the configurations for most of the temperatures studied. This can then be dealt 
with by importance sampling. However, with this algorithm, we encountered a nega- 
tive sign problem at the lowest temperature studied (/? = 96) for finite doping, and for 
higher temperatures, although the SDW persists for finite doping, we found no sign of 
development of a gap in the single-particle spectra, and no sign of the onset of or spin 
gap or building up of coherence in a staggered spin susceptibility, and no sign of the 
onset of superconductivity . This phase basically corresponds to an antiferromagnetic 
metal. In the DMFT based formalism, a metal-insulator transition is described as the 
onset of local or short-ranged gapless single-particle excitations and in this case, the 



transition is driven by carrier doping and accompanied by antiferromagnetic order. 1 27 



Since our primary goal is to investigate the relationship between short-ranged 
antiferromagnetic spin fiuctuations and superconductivity , we took a slightly different 
approach. The idea is to implement an algorithm such that once the system starts 
to build up an AF correlation (indicated by the appearance of a SDW as mentioned 
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above), two Ising variables are biased to stay at staggered configurations, namely 
(t, J,) or (i,t)- In other words, we attempt to build a machinery to do algorithmic 
projection onto low-lying states. One way to implement this idea is to use an identical 
random number array for both the A and B sublattices. More specifically, we first 
attempt a fiip of crA{Ti) and see if RA{Ti) > (random number) to accept the fiip and 
update. Subsequently, we do the same for the B-site using the same random number 
Xi- The motivation for this algorithm is basically a sublattice symmetry of our model 
which is manifested in the expressions for the accept ratio (please see the appendix 
for more details), i.e., 

RAin) = RB{Ti)\aB=-'TA,A^B,1^l (8) 

Once the system passes into an antiferromagnetically correlated critical domain (by 
which we mean a domain where the system acquires a tendency to pick staggered 
configurations) , the fiip which causes a transition from a parallel to a staggered con- 
figuration starts to get accepted more frequently than the fiip which causes the op- 
posite transition. For the sake of argument, let us consider the idealized case in 
which the former fiip is always accepted (that is to say, accept ratio for the former 
fiip is always greater than 1. Although this is generally not true due to the presence 
of fiuctuations, the argument below can be easily generalized to a more realistic case). 

Suppose the Ising variables for A and B sublattices are initially in the {a a, ctb) — 
(T) i) configuration at r — Ti. Then, due to the sublattice symmetry, if the fiip of a a 
is not accepted, the fiip of as is not likely to be accepted either, i.e., if i?A(Ti) < Xi 
for the fiip ((Ta,cb) : (T, i) ^ (i,i), then, it is likely that RsiTi) < Xi for the fiip 
((Ta, (Jfi) : (t, i) =^ (T, T) as well due to the sublattice symmetry. That is, if the initial 
configuration is a staggered one, a conditional probability that a ub fiip is not ac- 
cepted given that a a a is not accepted is high. This, of course, depends on the level of 
noise, thermal or quantum, already present, that is, the above statement holds in an 
averaged sense. The better a sublattice symmetry is preserved at each time slice, the 
stronger this correlation will be. Please note that if we use a random number j/j for 
the determination of a fiip of (Ts which is completely uncorrelated with Xj, then the 
probability of a ub fiip solely depends on the value of RB{Ti) and a rejection of a a fiip 
has no bearing on the ub fiip- If the fiip of a a is accepted and the system is updated 
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to (i, i) configuration, due to the presence of antiferromagnetic correlation, the fli p 
of ub will always be accepted(as assumed above). Therefore, in this case, whether or 
not to use the same random number will have no influence on the statistics. 

If the initial configuration is a parallel one, say, {(Ta,cfb) — (TjT)) then due to 
the presence of antiferromagnetic correlation, a a a flip will always be accepted and 
the system is updated to (i,t)- Therefore, in this case, it will make no difference 
whether or not to use the same random number for the (Tb flip- Thus, our algorithm 
is expected to bias the system towards staggered configurations more than a stan- 
dard algorithm, when a system acquires an antiferromagnetic correlation. The above 
argument can be easily generalized to the case where antiferromagnetic correlations 
are not ideal, i.e., the probability of a flip from a parallel to a staggered configuration 
is less than 1. 

Please note that, once the system starts to lose antiferromagnetic correlations 
(by which we mean that the transition probability between staggered and parallel 
configurations becomes equal), then the final population of the two configurations 
will converge to the same value. At this stage the effects of our algorithm will also 
start to diminish, simply because transitions from a staggered to parallel configuration 
become more frequent and the conditional probability discussed above becomes irrele- 
vant. Therefore our algorithm more or less approaches the standard Metropolis-style 
algorithm when the antiferromagnetic correlation diminishes, but starts to deviate 
from it when the system passes into an antiferromagnetically correlated phase. (Our 
algorithm respects spin rotational symmetry in the paramagnetic phase.) Please note 
that in our algorithm Ising variables can in principle take all possible configurations 
and no symmetry on the state is explicitly enforced. Therefore the symmetry and 
dynamics of the system still depends on the physics and remain nontrivial. 

As will be seen from results in Sec.III, our algorithm leads to the onset of coherent 
antiferromagnetic spin fluctuations for a certain parameter regime and unambiguously 
identifies a close tie between the onset of SC and of a coherent antiferromagnetic spin 
fluctuation. We will see that the results obtained by our algorithm, both for the 
dynamical and equilibrium properties of the model, are qualitatively consistent with 
the experimental findings for high-Tc cuprates. This suggests that the spin dynamics 
and its relation to the onset of superconductivity as found with our algorithm is more 
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or less consistent with what is reahzed in real systems. 



Once the imaginary-time correlation functions are evaluated, equilibrium quanti- 
ties can be directly calculated from the self-consistent Matsubara Green's functions, 
and dymanical quantities can be obtained from the via analytic continuation. We 
employed the Maximum Entropy(ME) ^ method to do this. For each set of 
parameters, we perform the iteration (3)-(5) and (7) until all the main features in the 
dynamical spectrum converge. In this study, we chose the number of time slices, L, 
to be such that At = ^ = 0.25 which should be sufficiently small to obtain a reliable 
value for the physical quantities of interest. [^. 



III. Results and Discussion 



A. Choice of Parameters 

We choose t = 0.5 throughout this paper such that the bare bandwidth of the under- 
lying Bethe lattice (equal to 2t) is the unit of the energy scale. U = 2.2 was chosen 
based on a preliminary study on the one-band Hubbard model where we found this 
makes the Neel temperature roughly optimal at half-filling. Then we sampled a few 
t' values to search for the one which leads to a Neel state with the experimentally 
observed < 5*2 >(~ 0.5/iB) at half filling. The optimal value for t' appears to be lo- 
cated somewhere between 0.1 and 0.2. In this paper, we show the results for t' = 0.2. 
(Setting t' = 0.1 does not modify the results in a qualitative way.) Although the 
above parameter setting has not been carefully adjusted to be the optimal choice, 
the overall agreement (qualitative and quantitative) with the experimental findings 
suggests that the above set of values are reasonable to refiect those for real materi- 
als. Again, let us emphasize here that our main interest is to study the qualitative 
relationship between superconductivity and antiferromagnetism. 
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B. Single-peirticle spectra 

We will present single-particle spectra both with and without the inter-cluster 
pairing coupling, and also compare them in the normal and superconducting phases. 
The single-particle spectra are defined as follows: 

G{uj) = -ImTrG{iuj) 

= G^L{^) + G^}A^) + G^U^) + G^U^) (9) 

Fig.2 shows the results at = 0.6(5 ~ 0.05) and = 0.3(5 ~ 0.17) without the inter- 
cluster pairing couphng. Superconductivity sets in at /3 ~ 16 for both dopings. The 
onset of superconductivity is determined by the onset of the average of the equal-time 
anomalous Green's function which is defined as: 

A single-particle gap starts to form when superconductivity sets in and a quasiparticle- 
like feature appears on top of the valence band as the temperature is lowered. Also, a 
dip which separates the quasiparticle peak and the relatively broad feature develops. 

Next, we study the case with inter-cluster coupling included where we have as- 
sumed gr'"*^" — g^nter (Jigcussed above, this heuristic extension of our 2-site cluster 
model is put in to give a more robust superconducting ground state, which allows us 
to extend our numerical calculations over a wider temperature range for studying the 
properties of the superconducting state. Fig. 3 shows a normal state single-particle 
spectrum for S ~ 0.05 and P — 16 which was obtained by setting all the anomalous 
components of the Green's function to zero (i.e., rj = 0). We confirmed that supercon- 
ductivity indeed develops at this parameter setting when the anomalous components 
are turned on. As we notice, in the absence of these components, the single-particle 
gap disappears and a sharp Kondo-like peak develops at the Fermi level, indicating 
that holes which participate in the pairing are strongly correlated. Fig.4 shows the 
the effect of inter-cluster pairing on the evolution of single-particle Green's function 
as a function of temperature for fj, — 0.7(5 ~ 0.04) and jj, = 0.3(5 ~ 0.17) where we 
have adopted the above values for t' and U for which g t as shown above. Although 
the overall features are similar to Fig.2, Tq is increased (superconductivity sets in at 
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/3 ~ 10) and the single-particle gap becomes more robust especially at smaller dop- 
ing. Experimentally, a single-particle gap has been found to decrease as a function of 
doping p2| . 



In all the cases shown above, superconductivity is found to set in before a robust 
gap has been formed (In particular, for /i = 0.3 with inter-cluster pairing). As we 
will show later, a signature of the onset of pairing fluctuations is more visible in the 
pair-pair correlation function. The onset of pairing fluctuations despite the absence 
of a robust gap is basically due to fluctuation effects which are neglected in the BCS 
mean-field theory. Here we are able to see this effect through our use of the dynamical 
mean-field theory approach. In order to confirm that the gap formation is due to the 
onset of superconductivity , we also studied the single-particle spectra in the normal 
state for the same parameter settings. Fig. 5 shows single-particle spectra at fixed 
doping {d ~ 0.04) as a function of temperature in the normal state. A "Kondo peak" 
is manifested in the spectra[^. For visual assistance, we overlay the results already 
shown in Fig. 5 (a) for the superconducting state. As the superconducting amplitude 
grows, the spectral weight near ep is depleted and transfered to the gap edges. These 
results unambiguously show that the single-particle gap is driven by the onset of su- 
perconductivity . Fig. 6 shows the single-particle gap (determined by the FMHW), Tc, 
and |F(0)p(please see Sec. 111(D)) as a function of doping at (3 = 16. Except for the 
low doping region, Ai particle roughly scales as which is what one would expect in 
the BCS picture. 

Next we compare our results with the ARPES data. Since the single-particle spec- 
tra calculated here are local, they should be in principle interpreted as angle-integrated 
quantities. Thus, only a qualitative comparison can be made with the angle resolved 
data. First, the quasiparticle peak as well as the dip for both doping is qualitatively 
consistent with the ARPES measurements. By comparing the spectra for the normal 
phase with those for the superconducting phase, the appearance of the quasiparticle 
peak in the superconducting phase can be qualitatively understood within the BCS 
picture in which quasiparticles in the vicinity of eF{\e — eF\ < Ai particle) in the normal 
phase form the Cooper pairs and the quasiparticle peaks appear at both gap edges. 
The appearance of the dip implies the separation of two energy scales, a sharp feature 
closer to ep and a broad feature below the sharp one. This broad feature is quite 
reminiscent of the "hump" which has been observed in the ARPES measurements. 
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Recently, the microscopic origin of the "hump" has been debated |^ Q , and argued 
that the underlying band structure of high-Tc cuperates is important for the onset 



of hump |33|. Here, we would like to make some remarks on this issue based on the 
results of our calculation. As shown above, the appearance of the hump is strongly 
correlated with the onset of superconductivity (or a quasiparticle peak). As we show 
later, a short-range antiferromagnetic spin fluctuation emerges as a coherent mode 
as superconductivity sets in. (Please note that our calculation is essentially exact 
and considers all energy scales involved). Therefore, it is reasonable to suspect that 
a hump-like feature is due to scattering of the photoelectrons by these coherent anti- 
ferromagnetic spin fluctuations. In fact, the separation of the quasiparticle peak and 
the hump-like feature is roughly proportional to the characteristic frequency of the 
antiferromagnetic spin fluctuations. This is an indication that the appearance of the 
two energy scales is due to a resonance. Please also note that since our calculation 
is fairly local, our results suggest that a peak-dip-hump feature seen in ARPES mea- 
surements does not depend on the details of the underlying band structure. 

The evolution of the single-particle gap as a function of doping is also consistent 
the ARPES experiments |^ ^ in which a pseudogap was found to monotonically 
decrease in the underdoped region and fall off more rapidly in the overdoped region 
as the doping increases. This suggest that Tc obtained in this study behaves more 
like T*. This observation is also consistent with the above finding that single-particle 
spectra start to show an anomaly at Tc. The actual ARPES spectra show a more 
pronounced quasiparticle peak than the one found here. This could be due to the fact 
that our DOS is an angle-integrated quantity and is also limited by the resolution of 
the Maximum Entropy algorithm. 



C. Two-Particle Correlation functions 

When a gap develops in a single-particle channel, a coherence or rigidity usu- 
ally start to build up in two-particle channels. A natural question is then if and how 
two-particle correlation functions show anomalous behavior as a result of the onset of 
superconductivity . This question is also motivated by the fact that the pseudogap 
has been observed in ARPES, ETM and STM (which basically measure single-particle 
spectra) and that INS, NMR and transport measurements (which basically probe the 
two-particle correlation functions in spin and charge sectors) have also found rich phe- 
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nomena as the pseudogap sets in. In this section, we present local susceptibilities of 
relevant bosonic degrees of freedom (both particle-hole and particle-particle channel) 
which appear to be playing an important role in high-Tc cuprates. 



1. Spin degrees of freedom 

Experimental studies of spin dynamics have been one of the most valuable source 
of information on the microscopies of high-Tc cuprates. In this section, we partic- 
ularly focus on resonance and a spin-gap features, and their relation to supercon- 
ductivity. The onset of the resonance which has been observed in inelastic neutron 
scattering (INS) measurements as a narrow peak at the antiferromagnetic wave vec- 
tor (5=(f , ^) (a: lattice parameter) and energy ~ 41mey(near the optimal doping) 
||37| , |38| , |40| , seems to be intimately related to the onset of superconductivity 
EO, B. The microscopic origin of this resonance as well as its connection to the onset 



of superconductivity has been recently debated by a number of authors JT^ |T8|, p3 
The spin gap is characterized by a depletion of low energy spectral weight in the 
dynamical spin susceptibility at commensurate wave vector Q and its onset appears 
to be correlated with that of a single-particle gap in ARPES spectrum The rel- 
evant quantity is the dynamic staggered spin susceptibihty Xz{Q,^) which basically 
contains the information about the collective excitations of spin degrees of freedom 
at wave vector Q. Using linear response theory, Xz{Q,^) is given as the following 
retarded staggered spin-spin correlation function: 



X.(g,a;) =-2 / dte''-'e{t) <[M,{t),MM]> 

J —oo 

where the staggered magnetization in our two-site cluster is defined as: 



Then, by the fluctuation-dissipation theorem and analytic continuation, Xz{Q,^) can 
be directly calculated from the imaginary-time staggered spin-spin correlation func- 
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tion which is given as: 

S,{Q,t-t') = < M,{t)M,{t') > 

Zmf 

The dynamic spin sucseptibihty, XziQ, can then be calculated using the fluctuation- 
dissipation theorem as follows: 

X.(Q,a;) = (l-e-^-)^,(Q,a;) 

Please see the appendix for more details of the calculation. 

One remark should be made here. Xz{Q,^) defined above should be interpreted 
as a measure of short-range fluctuations (i.e. the relevant length scale represented by 
the clustering) due to the two-site cluster nature of the model. That is, the above 
is a response function deflned for a probing fleld which acts only on the 
two sites in our cluster instead of the entire lattice. Thus, when Neel order sets in 
and the spin-spin correlation length becomes sufficiently longer than the cluster size, 
Xz{Q,^^) will start to deviate from the global staggered susceptibility defined as a 
response function of the entire lattice. Although we will not have precise information 
on the correlation length until we actually calculate it (probably by gradually increas- 
ing the cluster size which will become numerically expensive as mentioned above and 
is beyond the scope of this paper). Since a well-defined SDW is severely degraded 
as soon as the system is doped, and as mentioned earlier, the correlation length of 
spin fluctuations in high-Tc cuprates away from half-fllling as estimated by neutron 
scattering measurements is roughly of the order of a few lattice spacing, we believe 
that the local description given above should be a reasonable approximation. 

Fig. 7 and 8 show the cvohition of Xz{Q,^) as a function of temperature for 
different doping levels with and without inter-cluster pairing. In all cases, as the 
temperature is reduced, a feature starts to grow at some energy scale and becomes 
more coherent at lower temperatures. At the same time, low energy spectral weight is 
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gradually depleted and the spectrum eventually becomes gapped. These features are 
quite reminiscent of the observed behavior of the resonance and spin gap. Since the 
superconductivity starts to set in at /3 ~ 10 with inter-cluster coupling and at ~ 16 
without it, they appear to be correlated with the onset of superconductivity . The 
position of the resonance appears to be rather insensitive to temperature variation. 
Next, we plot XziQ, i^) as a function of doping at /3 = 16 for no inter-cluster coupling 
(Fig. 9 (a)) and with inter-cluster coupling(Fig9(b)). As long as superconductivity 
persists, a peak basically remains visible. The position of the peak appears to increase 
monotonically as a function of doping, and as a result, the spin gap also increases as 
the doping is increased (at least in the underdoped region). Upon passing into the 
overdoped region, the peak suddenly becomes incoherent and the spin gap disappears. 
We also found that where a resonance remains coherent, its energy is basically below a 
particle-hole continuum (i.e., u^es < 2Aj,c. See below for optical conductivity results). 
The overall qualitative feature of the evolution of Xz{Qi^) as a function of doping is 
basically similar for both cases. 



Now we compare the above results with the experiments. A resonance peak has 



been found by spin-flip INS in the superconducting phase [40, O, y]. The position of 



resonance Ures has been observed to soften as the doping is decreased from the optimal 
level [0, ^, ^ and appears to be relatively temperature insensitive |^, ^ ^ for 
a fixed doping. The magnitude of a spin gap also appears to decrease as the doping 
is reduced from optimal doping ^. These features are also qualitatively captured 
by our results. Furthermore, the relative energy scales among the resonance, spin 
gap and single-particle gap as found by INS and ARPES measurements are semi- 
quantitatively consistent with the results shown in Fig.4, Fig. 5, Fig.7 and Fig.8. 



To further investigate the relation of the spin spectra with superconductivity, 
we also studied the spectra in the normal phase. Fig.lO is XziQ,^^) for = 16 
and 6 ~ 0.17. For comparison, we also show the data in the superconducting phase 
with roughly the same parameter settings. In the normal phase, the coherence of 
the resonance is severely degraded and the existence of the peak can no longer be 
recognized. Also, XziQ,^^) acquires a substantial weight in the low energy sectors 
and the spin gap is completely filled up. These findings imply that the short-range 
antiferromagnetic spin fluctuations maintain their coherence by means of the pairing 
correlations. 
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2. Charge degrees of freedom 



Next, we present the dynamics of charge 



degrees of freedom. Since all the qualitative features found for the case of no inter- 
cluster pairing are also found when the inter-cluster pairing is switched on, we will 
only show the results for the latter case. Optical measurements have also revealed 
the emergence of a gap feature at some temperature T* which is higher than Tc in the 
underdoped region. A relevant physical quantity which probes the charge dynamics 
is the optical conductivity which is directly related to the retarded current-current 
correlation function. In our two-site cluster model, a local current operator can be 
defined in the following form: 



where the direction of current is that of intra-cluster hopping. This is basically a 
lattice version of current operator which takes a more familiar differential form in a 
continuum limit. Then, the Kubo formula for the paramagnetic part of optical conduc 
tivity aP^^"'{u) is given in terms of the following retarded current-current correlation 
function 



Since our clusters only have two sites, only the longitudinal part can be defined. 
As in the case of spin dynamics, we first calculate the following imaginary-time causal 
current-current correlation function: 



J 



a 




<j{'r)j{r') > 
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Ther detailed form of Xjj^r — r') is omitted here 0. Then, aP°-^°'{u) can be obtained 
by the fluctuation-dissipation theorem and analytic continuation: 



We plot /2e[cTP"™(cj)] for /i = 0.3 and /i = 0.7 in Fig.ll. For fx = 0.7(Fig.ll(a)), 
at higher temperatures, a Drude-like peak is visible due to the fact that the single- 
particle spectrum has a nonzero weight at ei? (Please see Fig.4(a)). As the temperature 
is lowered, the Drude-like weight diminishes and eventually a gap opens up as super- 
conductivity develops. Even when the system passes into the superconducting phase, 
since the single particle spectrum does not become fully gapped, a small weight still 
remains in the low energy sectors above P = 16. Again, this is due to fluctuation 
effects which are not included in the BCS theory. The size of the gap is quantitatively 
consistent with a single-particle gap (i.e., AopUcai — '^'^i particle)- At this doping, a 
charge gap feature (i.e., transition between the lower and upper Hubbard bands as 
seen in Fig. 11 (a)) is visible for all temperatures. A feature corresponding to the quasi- 
particle peaks in the SC state becomes visible as the temperature is lowered(/3 = 24). 
In Fig.ll(b)(for fx = 0.3), a large Drude-like peak appears at high temperatures due 
to a substantial weight at ep in the single-particle spectrum(Fig.4(b)). As the tem- 
perature decreases, the weight in the low energy sectors disappears and a gap feature 
develops. Similarly, due to residual single-particle weight at ep, a robust gap does 
not develop in charge spectra until the temperature is lowered to P = 24. Please note 
that a feature for the charge gap for fi = 0.3 is hardly visible, which is consistent 
with our single-particle spectra at this doping(Fig.4(b)). This is basically due to a 
low electron density. An addtional feature on the quasiparticle peak (i.e., the feature 
at the gap edge), however, is much more robust at /i = 0.3 (as can be expected from 
Fig.4(b)) and starts to show up already at /3 = 16. This feature on the gap edge has 
been observed in the infrared spectra for the optimally doped both single and bi-layer 
high -Tc materials. The size of the gap for fi = 0.3 is reduced from that for fx = 0.7 
and is also quantitatively consistent with the single-particle gap as shown in Fig.4(b). 
Please note the wide range of energy scales over which the weight is redistributed 
as a function of temperature. This strongly temperature-dependent spectral weight 
has been observed in the optical conductivity spectrum of high-Tc cuperates |Q and 
suggests that the charges participating in the pairing are strongly correlated. 
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What we showed (numerically) in the above is basically that a paramagnetic cur- 
rent disappears in the superconducting phase below a certain energy scale(~ ^optical)- 
Therefore, this will naturally lead to the onset of the Meissner effect caused by a resid- 
ual diamagnetic response which comes from the mobile Cooper pairs. The magnitude 
of the diamagnetic response is proportional to the Cooper pair density. 

3. Pair degrees of freedom We now turn to the particle-particle channel. 
The quantity of interest is the pair-pair correlation function. This function basically 
contains information on the amplitude fluctuations of the superconducting order pa- 
rameter. A local d-wave - like pair-pair correlation function Xd{'^) can be defined in 
the following way: 



where a local d-wave order parameter centered at a particular A-site is defined as: 



Please note that each pairing for a given nearest-neighbours appears as a singlet which 
is a direct consequence of taking a zero center-of-mass momentum of the pair (please 
see below) . This expression can be obtained by simply transforming back to real space 
the following more familiar definition of d-wave operator in 2d momentum space: 




- t^Ai(^)4iTW-^\W42TW 

+ Al{t)^BZ^{t) + ^Al{t)^BA^{t) 



{C0s{p:c) - COS 



{Py))dl^{t)dl,^l{t) 



P=(Px,Py) 



24 



After performing an inverse Fourier transformation, this leads to the following (global) 
d-wave - like operator which is basically a sum of the above local d-wave operator 
over the all A-sites: 



The choice of A-site as a reference is arbitrary and B-sites could be chosen just as 
well. 

Since the clusters only have two sites, we need to make an assumption about the 
symmetry of the dynamics, namely a rotational symmetry among the four nearest- 
neighborhoods. The above Xd{^) is defined with respect to a particular A-site and 
contains the pairing correlation between its four nearest-neighboring B-sites. Here, 
we basically make the assumption that the dynamics of these four pairings are the 
same except for their phases. We view this as a way to represent the effects of inter- 



cluster correlations within the constraints of our model [0. Therefore, the above 
local Xd{A,u!) can be collapsed into a single two-site cluster. Then, we calculate its 
corresponding imaginary-time causal correlation function: 
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+ 5rr' < G'2f\r, r) + G^^'V, r)) >.(o (10) 



where Xdi^ — t') is an average per pair. Basically, the product of F's and G's cor- 
respond to a (local) pair hopping and pair breaking fluctuations, respectively. After 
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the analytic continuation (Maximum Entropy), we finally obtain a local dynamical 
pair-pair correlation function Xd{^^)- 

We plot Xd(^) for fi = 0.7(Fig.l2(a)) and /i = 0.3(Fig.l2(b)). We see that for 
both dopings, as the temperature is lowered, a narrow feature grows around a; = 
which indicates that the superconducting order parameter develops a well-defined 
static component. The temperature at which this narrow feature develops basically 
corresponds to determined from the onset of instantaneous nearest-neighbor pair- 
ing correlation. 

From a phenomenological point of view, pair hopping and pair breaking fiuctua- 
tions correspond to short-wavelength phase and amplitude fiuctuation of the super- 
conducting order parameter, respectively. The latter is more intuitively transparent 
but the former becomes readily clear if one notes that: 



Therefore, we can see in more detail the nature of the onset of a superconducting 
phase by studying both pair hopping {F x F terms in (15)) and pair breaking {G x G 
terms in (15)) fiuctuations. From detailed studies we found that pair breaking 
fiuctuations are more strongly gapped at Tc for /i = 0.7 than for // = 0.3, and this 
implies that the amplitude fiuctuations are more enhanced for /i = 0.3 when super- 
conductivity switches on. These findings basically suggest that Cooper pairs are more 
tightly bound in the low doping region (which is consistent with the doping evolu- 
tion of the single-particle gap). Since phase fiuctuations are typically energetically 
cheaper than an amplitude fiuctuations (this will typically be a gapless Goldstone 
mode in the long wavelength limit), a remedy for the over-estimate of Tc in the low 
doping region may be sought by incorporating a phase fiuctuation mechanism, or 
equivalently charge localization in the dual picture. (In the present model, the coher- 
ent pair-hopping fiuctuations are gapless as seen in Fig. 12, and phase coherence can 
not be suppressed, i.e., a charge localization mechanism is absent.) Although such 



{pair hopping) 




where 




(11) 
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localization mechanisms have been proposed |i6, |^|, phenomenology in the lightly 
doped region is still under debate Please note that the physics which governs 
global phase coherence is a low energy long-wavelength phenomenon (depending on 
the localization length of Cooper pairs) and is thus not possible to address by the 
present approach based on the small cluster size. 



D. Equillibrium properties and phase diagram 



Since our main focus is antiferromagnetism and superconductivity , we inves- 
tigate those two phases. We show our results both with and without the inter-cluster 
coupling. 

An antiferromagnetic phase is signaled by the onset of the thermal average of a 
staggered magnetization induced by an infinitesimal staggered field. Here, instead, we 
chose to induce symmetry breaking by appropriately initializing the Ising variables 
for each sublattice. Although the staggered magnetization is in principle a vector 
order parameter a natural choice is its z-component, since the Hubbard- Stratonovich 
decomposition was performed along the z-axis. This is given as: 



<S\> 



<S\> 

< flA^ > 

< flB^ > 



<s%> 

< flAi > 
<hBl> 



where < ... > means a thermal average. In Fig. 13 (a), the staggered magnetization, 
< >, is plotted at half-filling as a function of temperature. Since superconductiv- 
ity does not set in exactly at half filling, there is no distinction between the two cases. 
The staggered magnetization disappears roughly at = 6. Upon doping, < > 
almost suddenly drops and disappears before S reaches ~ 0.009. A detailed study of 
the destruction of antiferromagnetism as a function of doping is outside the scope of 
this paper, and we did not sample sufficiently many doping levels to determine the 
nature of transition, such as the order of transition. 

Similarly, if we employ the definition of the BCS picture, an onset of supercon- 
ductivity can be defined as the onset of an average of equal-time anomalous Green's 
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function: 



All the on-site pairing components, F\\^{0), ^^^(O), F^^(O) -^_bb(0)5 always 
one to three orders of magnitude smaller. We turn on a small 1] and see if the system 
can sustain a finite |F(0) p. Therefore, according to linear response theory, the appear- 
ance of a finite |F(0)p basically corresponds to divergence of the pairing fluctuations. 
In Fig. 13, |F(0)p is also plotted as a function of doping for /3 = 16(Fig.l3(b)) and 
as a function of temperature for large doping(Fig. 13(c)) and small doping(Fig. 13(d)). 
We found that the inclusion of inter-cluster coupling systematically enhances Tc|p. 
The doping which gives a maximum |F(0)p seems to be somewhere between /z = 0.4 
{6 ~ 0.13) and fi = 0.3(5 ~ 0.17) for both cases. A determination of the doping 
which gives a maximum Tc requires a detailed study involving varying both temper- 
ature and doping and was difficult to identify with the same accuracy(it appears to 
be located at least between 6 ~ 0.009 and 6 ~ 0.17). Note that a superconducting 
phase appears as soon as the system is doped. 

In Fig. 14, we show the phase diagram as a function of doping and temperature in 
the presence of inter-cluster coupling. (A phase diagram without inter-cluster pairing 
is also shown here. The overall topography is quite similar Q.) Note that the overall 
topography is quite consistent with that of high-Tc cuprates: a sudden suppression 
of AF order away from half-filling, an upper critical doping for superconductivity at 
roughly 6c ~ 0.3 and an optimal doping is located 0.13 < 6 < 0.17. For t' = 0.1, all 
the qualitative features are essentially the same although the upper critical doping 
and the staggered magnetization at half-filling appear to increase slightly. 

One notices, however, that a superconducting amplitude is quite alive even in a 
lightly doped region, i.e., the lower critical doping is essentially zero. This feature 
was found for a related model in a BCS mean-field approximation (which can 
be thought of as a static version of the present approach) and later interpreted as 
the onset of a phase stiffness(i.e., the amplitude part) of the superconducting order 
parameter |^ which exactly corresponds to |F(0)p here. A locally constructed su- 
perconductivity was also found to coexist with antiferromagnetism for a relatively 
wide range of doping in the dynamical mean-field theory study of equilibrium state 
of a cluster model in which local dynamics are more constrained [jTB[. As mentioned 
in the previous section, the physics of underdoped region is still under debate. 
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IV. Summciry and Conclusions. 



As mentioned previously, the relatively short coherence and spin-spin correlation 
lengths which are observed in the in high-Tc cuprates constitute the key motivation 
of the present study as well as underpinning the basis of our approach. The subject 
of this paper is to investigate the role of short-range quantum fluctuations in high-Tc 
superconductivity and, in particular, the interplay between pairing and spin fluctu- 
ations. In order to achieve a more complete understanding of the strongly coupled 
dynamics whose calculation goes beyond the bounds of BCS mean-fleld we apply dy- 
namical mean-fleld theory to a simple 2-atom cluster model. We showed that this 
model, despite its simplicity, not only reproduces the main characteristic features 
of high-Tc cuprates(i.e., the basic topography of the phase diagram, a quasiparti- 
cle feature in single-particle spectra, the spin-gap, neutron resonance, d-wave- like 
pairing, etc) with qualitative and semi-quantitative consistency, but also revealed a 
strong correlation between coherent short-ranged antiferromagnetic spin fluctuations 
and pairing correlations: i.e., if superconductivity is suppressed the antiferromagnetic 
spin fluctuations lose their coherence. Thus our approach allows us to directly study 
the relation between the dynamics of the spin and pair degrees of freedom. 

We would like to emphasize that the present model was constructed from strongly 
correlated electron degrees of freedom only and no prior assumption was made for the 
existence of any intermediate energy scale which typically represents some sort of com- 
posite bosonic degrees of freedom in either particle-hole or particle-particle channel. 
In the present formalism, bosonic degrees of freedom arise naturally as a consequence 
of gap formation in the single-particle sector. In our model, away from half filling, 
coherence builds up in a singlet particle-particle channel and this is accompanied by a 
development of coherence in the staggered spin channel at some characteristic energy 
scale. The only nontrivial assumption that we put into the present model is basically 
that the correlation length of spin and pairing fiuctuations is sufficiently short and 
that antiferromagnetic fiuctuations are a relevant spin fiuctuation mode to consider. 
In this context we have been able to show that the spin fiuctuation resonance which 
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has been discussed by a number of authors [17, as a spin - flip exciton-hke mode 
is a generic feature of our model. This suggests that it survives beyond the RPA 
approach and should not depend on details of the band structure of the system. Thus 
it appears to be an intrinsic property of the strongly coupled superconducting state 
in the presence of strong on-site repulsion. 



From our results on the single particle excitations of the model, and the fact 
that we can correlate the "hump feature" seen in ARPES spectra with that seen in 
our model, we conclude that the spectral weight of the spin fluctuation resonance is 
sufficient to show up as an energy loss peak in the photeoemmission spectra. This is 



in contrast to recent argument sp4|, |33[| suggesting that details of the band structure 
are important in this interpretation of the ARPES data. From the point of view of 
the local physics represented in our model, it appears that this feature is an intrinsic 
property of Hubbard-type models of superconductivity. 

The present approach is essentially an analogue of Weiss mean-field theory in 
which fiuctuations along the imaginary time are exactly treated but all the spatial 
fiuctuations whose wavelengths are longer than the size of the cluster are averaged 
out. Unlike theories which are based on momentum space representation, it is based 
on a local picture. Therefore, when the physics in the thermodynamic limit is domi- 
nated by low-lying long-wavelength modes, this local description is expected to break 
down. (When the correlation length becomes sufficiently longer than the cluster size 
but still finite, the present model will incorrectly assume that true long-range order 
is established). Typically, it is the dimensionality of the system which plays a crucial 
role in deciding the fate of the thermodynamic limit when long- wavelength fiuctua- 
tions dominate. (Of course, the length scale of interactions among the local order 
parameters, short-ranged or long-ranged, is also important and this is assumed to 
be of short-range in the present study). Therefore, physics which depends on the 
dimensionality in an essential way will be difficult to address in this approach. (This 
approach is essentially a "0 + 1" dimensional formalism.) Instead, the present model 
is looking at the portion of physics which is rather insensitive to the dimensionality, 
i.e. the short-range short-wavelength modes. Such features are relatively dimension 
independent - indeed, superconductivity and similar anomalous features have been 



reported in quasi- ID doped Hubbard ladder systems as well [H9 
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Before we conclude, we note that recent studies [13, IG] have looked at dynamical 
properties of four-site cluster models. Although these authors have not yet investi- 
gated the excitation spectra of such models, we expect them to lead to qualitatively 
the same conclusions that we have reached based on our 2-site cluster model. A four- 
site square cluster model would also allow consideration of other kinds of long range 
order such as a chiral flux phase A four-site cluster model would also allow 

for study of the dynamics of 7r-operators in the presence of broken symmetry phases 
(effectively in a thermodynamic limit) and the connection between SO (5) symmetry 
fS^ , ^ and superconductivity . This would be complimentary to results of exact 
numerical workp4| for finite clusters. 



The present approach is to transform the original lattice problem to an effective 
self-consistent local which can be exactly solved in principle. Thus, it can provide 
essentially exact solutions to the original problem within the range of validity of this 
transformation. Certain aspects of high-Tc superconductivity appear to fit reason- 
ably well into the regime in which this transformation is valid. The present model is 
perhaps the simplest nontrivial one in the context of a self-consistent cluster model in 
the dynamical mean-field theory approach. Yet, it appears to succeed in addressing 
some nontrivial aspects of strongly correlated electron systems which would be diffi- 
cult to approach by other means. The results we have presented clearly suggest that 
as long as the physics at hand is of short-range this approach can be quite effective 
and promising. 



In summary, we constructed a self-consistent two-site cluster model in the dy- 
namical mean-field limit in which short-range short- wavelength fiuctuations of charge 
and spin degrees of freedom are treated exactly in the presence of superconductivity. 
The equilibrium properties of our model (superconductivity and antiferromagnetism) 
as a function of doping and temperature reproduce the overall qualitative and ba- 
sic quantitative features of the phase diagram of high-Tc cuprates. The behavior of 
single-particle and two-particle spectra of our model in the superconducting phase 
can be interpreted to give a consistent account of anomalous features of high-Tc 
cuprates such as the pseudogap (charge and spin) and resonance observed in ARPES, 
INS, NMR and other optical measurements. 
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APPENDIX A: Complete self-consistency equation 

If we systematically integrate out all the hgand degrees of freedom, we obtain 
the following self-consistency condition which constitutes the lattice nature of the 
problem. The physical meaning of each term can be diagramatically illustrated in t 
he same manner as we showed in sec.III: 



QBsi^^n) = il^n - ed + fJ- - - — - fG^2A{'i'^n) 

GaA (i^n) = iuJn + ed- f^ + K + — - t^G^^B (i^n) 

~ ^AB^BA^BBy-^n) 

Gbb {i^n) = ii^n + ^d- 1^- K + ^ - fG^J:j^{iuJn) 

~ ^ ^BA'iAB'^AAV^n) 



GABi^^n) = t' -Mj^^izUr^) 

~ ^^iAB^AB^BAi^'^n) 

+ t\(^}BF^Jj,{iu;n) + i^2BFl\{i^n)) 
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OBAii^n) 



= t' - t'^G^^Bii^n) 

~ t'^^BA^^BA ^AB {^^n ) 



<,BA^ AB 



= t'-ec 



u 

BA 



.lOJr. 



~ ^'^i^AB^AB G^BA (^"^n ) 
~ ^ 'iBA'iBA^ ABV'^n) 

- {Cba Fab (^^n ) + Cba ^ab ' 



iiOn)) 



iUn)) 



^AAii^n) = t^Fj^Bii'^n) 

- t'Ci^B^WBii^n) 

- t'{eAWUi^n)-dWJB{iu^n)) 

- t\ejAG^^Bi^^n) - ^^JbGU^^u)) 

PJeii^n) = t^F\i^{iUn) 

- t'{dWL{i^n)-eA^BG^2A{i^n)) 



- t'{C\'BrFi\{zU^n) 

- t'^!^BiG\^Ai^^n) - G^U^^^)) 
^BAii^n) = -V + t^F^liioOn) 

- t'iejAfFUiUr.) 

- t'^^JAiG^Bii^n) - G^siicOn)) 
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- t'iC'LfFl'Ai^U^n) 



APPENDIX B: Manifestation of a sublattice symmetry in the formula 
of the accept ratio 

Since the Nambu-Gorkov representation is used in (9), we first need to trans- 
form to tlie usual representation. This can be done by transforming only the down 
spin components in the following way: 



This will transform the accept ratio to the following form: 



GAA{n,ri) I — > 1 - G^J;j^{Ti,Ti 
G^jB{n,ri) I — ^ 1 - G\^B{n,ri 



XaAin) _ p-Ao-A(-ri)\2 pT-L 
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in which the sublattice symmetry is manifestly apparent. 
APPENDIX C: Calculation of x%{^) 

In this appendix, we show the details of the procedure to compute a staggered spin 
susceptibility X(^{^) from an imaginary-time causal staggered spin-spin correlation 
function which is given as: 

S.{Q,t-t') = < M,{t)M,{t') > 

Zmf 

The above trace involves four point correlators and a Wick contraction needs to be 
performed. Since the effective action is bilinearized by the Hubbard- Stratonovich 
transformation, the functional integral over the fermionic degrees of freedom can be 
exactly performed and the above trace simply becomes a sum of the products of local 
Green's functions averaged over the Ising variable. In other words, all the vertex 
corrections at the local level are decoupled by means of Hubbard- Stratonovich trans- 
formation and re-absorbed into the dependencies of the local propagators on the Ising 
variable. For example, the term < nAi{T)nBi{T) > can be integrated in the following 
way: 



-kG^b {'r,T)>a{i) 

where < ... >a{i) means an averaging over the Ising variable. The normal order for 
down spin components is reversed due to the Nambu representation. The full expres- 
sion of Sz{Q,t) becomes quite lengthy and is omitted here 0. The averaging over 
the Ising variable is performed by the QMC sampling. Since Sz{Q,t) is a scattering 
function(i.e., fluctuation of a spin degrees of freedom), the transfer function which 
defines the relation between the imaginary-time and real-time quantities becomes: 
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1 /• + 00 

S^{Q,r) = — due~^'^S,{Q,uj) 

ZTT J-oo 

where Sz{Q, to) must satisfy the equation of detailed balance Sz{Q, —00) = e'^'^SziQ, uj) 
(Please note the difference from the fermion operators (10).) Here, the equation of 
detailed balance can be used as a constraint for the default model in Maximum En- 
tropy or simply absorb it into the transfer function, which is the approach we 
adopted here. 
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Figure 1 : A schematic picture of our model. Elliptic curves define the two-site 
cluster. All sites have the same nonzero U, all inter-cluster hopping have the same t, 
and all intra-cluster hopping have the same t'. 

Figure 2 : Evolution of single-particle spectra as a function of temperature 
without inter-cluster anomalous terms, (a) is for a small doping(5 ~ 0.05) and (b) 
for a large doping(5 ~ 0.17). 

Figure 3 : Single-particle spectrum in the normal phase for P — 16 and 5 ~ 0.05. 

Figure 4 : Evolution of single-particle spectra as a function of temperature with 
inter-cluster anomalous terms, (a) is for small doping(5 ~ 0.04) and (b) for large 
doping (5 ~ 0.17). 

Figure 5 : DOS of normal(solid line) and SC(dashed hne) state for /x = 0.7 as 
a function of temperature. A gap feature clearly builds up as the superconducting 
correlation develops. 

Figure 6 : A single-particle gap(diamond), Cooper pair density(oc |F(0)p, cross) 
and Tc(square) vs. doping at P — 16. Vertical units are arbitrary. 

Figure 7 : XziQ,^^) as a function of temperature for the underdoped region. 
(a):without inter-cluster pairing and (b):with inter-cluster pairing. 

Figure 8 : XziQ,^^) as a function of temperature near the optimal doping region. 
(a):without inter-cluster pairing and (b):with inter-cluster pairing. 

Figure 9 : Xz{Q:^) as a function of doping. (a):without inter-cluster pairing 
and (b):with inter-cluster pairing. 

Figure 10 : Xz{Q:^) spectra of normal(dash line) and SC(solid hne) for 13—16 
and II — 0.3. Coherence is severely degraded and the gap is filled in the normal phase. 

Figure 11 : i?e[cr^"''"(a;)] as a function of temperature for /j, — 0.7(a) and fx — 
0.3(b). The horizontal unit is 2t and the vertical unit is arbitrary. 
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Figure 12 : Xd{^) as a function of temperature for fj, — 0.7(a) and fx — 0.3(b). 
The horizontal unit is 2t and the vertical unit is arbitrary. 

Figure 13 : (a)M2 vs. temperature at half filling. Units are fXB and 2t for 
vertical and horizontal axis, respectively. (b)|F(0)p vs. doping at /3 — 16. A unit for 
the vertical axis is arbitrary. |-F(0)p vs. temperature for fj, — 0.7(c) and fj, — 0.3(d). 

Figure 14 : Schematic phase diagram obtained based on our model. The onset 
of SC was determined by the appearance of a thermal average of an instantaneous 
pairing amplitude. Error bars correspond to the size of increments of sampling points 
we studied. 
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Figure 1 : A schematic picture of our model. Elliptic curves define the two-site 
cluster. All sites have the same nonzero [/, all inter-cluster hopping have the same t, 
and all intra-cluster hopping have the same t'. 



42 



(a) \i = 0.6 

1 1 I 

p = 24 




-1.5 -1 -0.5 0.5 1 1.5 



0) 



(b) [I = 0.3 




-1.5 -1 -0.5 0.5 1 1.5 



CO 

Figure 2 : Evolution of single-particle spectra as a function of temperature without 
inter-cluster anomalous terms, (a) is for a small doping(5 ~ 0.05) and (b) for a large 
doping((5 ~ 0.17). 
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Figure 3 : Single-particle spectrum in the normal phase for j3 — 16 and S ~ 0.05. 
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(a) n = 0.7 




(b) n = 0.3 
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Figure 4 : Evolution of single-particle spectra as a function of temperature with 
inter-cluster anomalous terms, (a) is for small doping(5 ~ 0.04) and (b) for large 
doping ((5 ~ 0.17). 
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Figure 5 : DOS of normal(solid line) and SC(dashed line) state for // = 0.7 as a 
function of temperature. A gap feature clearly builds up as the superconducting 
correlation develops. 
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Figure 6 : A single-particle gap(diamond), Cooper pair density(oc |F(0)p, cross) 
and Tc(square) vs. doping at P — 16. Vertical units are arbitrary. 
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Figure 7 : XziQ,'^) as a function of temperature for tlie underdoped region. 
(a):without inter-cluster pairing and (b):with inter-cluster pairing. 
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(a) No inter-cluster pairing (|a = 0.3) 




(b) Inter-cluster pairing ^ = 0.3) 




Figure 8 : XziQ,'^) as a function of temperature near optimal doping region. 
(a):without inter-cluster pairing and (b):witli inter-cluster pairing. 
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(a) No inter-cluster pairing (P = 16) 




(b) Inter-cluster pairing (P = 16) 
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Figure 9 : XziQji^) as a function of doping. (a):without inter-cluster pairing and 
(b):with inter-cluster pairing. 
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P = 16, 5-0.17 




Figure 10 : XziQi'^) spectra of normal(dash line) and SC(solid line) for (3 = IQ and 
= 0.3. Coherence is severely degraded and the gap is filled in the normal phase. 
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Figure 11 : RelaP""^"" (u)] as a function of temperature for /j, 
The horizontal unit is 2t and the vertical unit is arbitrary. 



0.7(a) and ii = 0.3(b). 
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Figure 12 : Xd(a;) as a function of temperature for fi = 0.7(a) and fi = 0.3(b). Tlie 
horizontal unit is 2t and the vertical unit is arbitrary. 
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(a) 6 = (b)P = 16 




Figure 13 : (a)M2 vs. temperature at half filling. Units are hb and 2t for vertical 
and horizontal axis, respectively. (b)|F(0)p vs. doping at /? = 16. A unit for the 
vertical axis is arbitrary. |F(0)p vs. temperature for ji = 0.7(c) and ji = 0.3(d). 
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Figure 14 :Schematic phase diagram obtained based on our model. Tlie onset of SC 
was determined by tlie appearance of a tliermal average of an instantaneous pairing 
amplitude. Error bars correspond to the size of increments of sampling points we 
studied. 
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